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Abstract: We explain how angular correlations in leptonic decays of vector bosons 
and top quarks can be included in Monte Carlo parton showers, in particular those 
matched to NLO QCD computations. We consider the production of n pairs of 
leptons, originating from the decays of n electroweak vector bosons or of n top quarks, 
in the narrow-width approximation. In the latter case, the information on the n b 
quarks emerging from the decays is also retained. We give results of implementing 
this procedure in MC@NLO. 
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1. Introduction 

Accurate predictions for the spectra of the leptons emerging from decays of vector 
bosons or top quarks are important for a variety of studies at hadron colliders, such 
as acceptance computations, tests of QCD, and searches for new physics. Theoretical 
computations should be based on all Feynman diagrams in which the corresponding 
leptons are external legs. In general, not all such diagrams are resonant diagrams, i.e. 
those in which the leptons directly emerge from a vector boson propagator (which, 
in the case of top decays, in turn is directly connected to the top quark via a Wtb 
vertex). Usually, however, predictions based on computations that retain only the 
resonant diagrams are excellent approximations to those based on the fuller set of 
diagrams, owing to the rather narrow widths of the vector bosons and top quarks; 
the more so in the presence of final-state cuts which are designed to enhance on-shell 
contributions. 

A further approximation can be made, which we call the decay chain approxi- 
mation: resonant diagrams are replaced by diagrams relevant to the production of 
on-shell vector bosons or top quarks, times the diagrams corresponding to the matrix 
elements for the decays. In this way, off-shell effects are lost, but they can be recov- 
ered to some accuracy by reweighting the results of the decay chain approximation 
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by a Breit-Wigner function. There is another piece of information that is lost in 
the decay chain approximation, and cannot be recovered, namely that on production 
angular correlations (more precisely, angular correlations due to production spin cor- 
relations). Let us denote by P the decaying particle (a vector boson or a top in our 
case), and by di, . . ., d n its decay products, and consider the hard process 

a + b — > P(^rf! + --- + rf n )+X, (1.1) 

with X a set of final-state particles which may also contain other decaying vector 
bosons or top quarks. The process of eq. (1.1) is said to have decay angular correla- 
tions if the matrix elements of the corresponding resonant Feynman diagrams have 
a non-trivial dependence 1 on (di-dj). Clearly, decay correlations are always present 
if the particle P has spin different from zero. The process of eq. (1.1) has production 
angular correlations if its matrix elements have a non-trivial dependence on (d{-a), 
(di-b), or (dfX). It is therefore clear that the decay chain approximation can account 
for the decay correlations, but not for the production correlations. 

The decay chain approximation has obvious advantages, leading to much simpler 
computations (especially at higher orders) owing to the reduced multiplicity of the 
final state. Still, it is not acceptable if the spectra of the decay products must be pre- 
dicted with some accuracy. The aim of this paper is to introduce an approach to the 
computations of lepton spectra as given by resonant diagrams, which uses the decay 
chain approximation but also correctly accounts for production angular correlations. 
The method is primarily intended to be applied to parton shower Monte Carlos, 
including those that implement NLO QCD corrections such as MC@NLO [1, 2] or 
POWHEG [3]. The idea stems from the following observation: the matrix elements 
computed with the resonant diagrams are bounded from above by the matrix ele- 
ments obtained by eliminating the decay products and putting the parent particles 
(vector bosons and/or top quarks) on-shell, times a process-independent constant. 
One can therefore use the latter matrix elements (which we call undecayed matrix 
elements) to perform computing-intensive tasks for which production correlations 
are not an issue. When the four-momenta of the parent particles are available, the 
resonant diagrams (we refer to the corresponding matrix elements as leptonic ones) 
are used in the context of a simple hit-and-miss procedure to generate the leptonic 
four-momenta. 

In order to apply a hit-and-miss procedure, we need upper bounds on the decay 
matrix elements that are universal with respect to the production process. These 
are derived in the following section, first for vector boson, then for top quark de- 
cay, and finally for final states containing several vector bosons and/or top quarks. 
The practical application of these results is discussed in section 3. The inclusion of 
angular correlations in NLO computations is hampered by the presence of virtual 

1 We denote here a particle and its four- momentum by the same symbol. 
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corrections and the necessity for subtraction terms, which mean that one has to deal 
with expressions that are not simply matrix elements squared, and therefore are not 
necessarily positive-definite. This implies that the scheme we propose in this paper 
is such that angular correlations are not accurate to NLO in the whole phase space, 
but are correct to NLO for hard real emissions and to LO in soft and collinear re- 
gions. Obviously, one can implement angular correlations exactly to NLO accuracy 
by using lepton matrix elements in all the steps of the computation. In this paper, 
however, we are solely interested in the decay chain approximation. Illustrative re- 
sults of our approach, obtained with MC@NLO, are presented in section 4, followed 
by our conclusions in section 5. An appendix presents an alternative derivation of 
the upper bound for vector boson decay, which may clarify some of the assumptions 
involved. 



2. Upper bounds for the leptonic matrix elements 

In this section, we derive the universal factors that, when multiplied by the unde- 
cayed matrix elements, give an upper bound for the leptonic matrix elements. After 
introducing some notation, we shall treat the cases of the vector bosons and of the 
top quarks in turn. 

2.1 Notations 

We shall always denote by 

V — > 11 (2.1) 

the decay of the vector boson V = W or Z into a lepton-antilepton pair, which 
means that in the case of W decay 1 is not the antiparticle of /. In our conventions, 
the Vli vertex is 

-iF vl » (V vl - A vn5 ) , (2.2) 



where 

9w 



2 cos 9 W ' 



Vzi = V h A Z i = Au (2.3) 



F w = ^=, V wl = l, A wl = l. (2.4) 
2y2 

We shall consider the process 

a(Pi) + b(P 2 ) — > Vi (gi) + . . . + V n (q n ) + X(x) (2.5) 

— > k(h) + h{k 2 ) + ... + Z n (A; 2n _i) +J n (k 2 n) + X(x), (2.6) 

where 

Qi = hi-i + hi , (2.7) 
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and X collectively denotes any particles not originating from a vector boson decay. 
It is particularly convenient to write the phase space of the final-state particles of 
eq. (2.6) as follows 

n dq 2 
d<S> n+1 *(P 1 + P 2 ;q 1 ,...,q n ,x) d<S> 2 (qf, hi-i, hi) ■ (2.8) 

i=l 

As the notation 1* suggests, we treat the particles X as a single particle with mass- 
squared x 2 and four-momentum x, since the individual four-momenta of the particles 
X are irrelevant in what follows. On the r.h.s. of eq. (2.8), the two-body phase spaces 
account for the decays 

VMi) — > h(k2i-i) + h(k2i). (2.9) 

The factorization formula of eq. (2.8) is exact: the vector bosons are off-shell, and 
their virtualities q 2 (i.e., the invariant masses of the lepton pairs) are explicitly in- 
tegrated over. This decomposition has an obvious physical interpretation in the 
context of resonant diagrams. 

In the case of processes involving top quarks, we shall deal with 

a(P0 + b(P 2 ) — > tM + ... + t n ( Pn ) + X(x) (2.10) 
— > ^i(gi) + 6i(ri) + • • • + W n (q n ) + b n (r n ) + X(x) (2.11) 
— > h(ki) + vi(k 2 ) + foi(ri) + . . . + l n (k 2 n-i) + v n {k 2n ) + b n (r n ) 
+X(x) , (2.12) 

where t can be either a top or an antitop. As in eq. (2.8), we can also write the exact 
phase-space factorization 

d$ 3n+ i*(Pi + P 2 ; h,..., k 2n ,ri, ...,r n ,x) = 

d$ n+ i*(Pi + P 2 ;pi, . . • ,p n ,x) [[d^3(pi;k 2i - 1 ,k 2i ,r i ) , (2.13) 

i=l 

with the three-body phase spaces on the r.h.s. accounting for the decays 

kipi) — > WM + biin) — > liifai-J+Viifad + biiri). (2.14) 
2.2 Vector boson decay 

We start by considering the production of one U pair, and we neglect the Zj 7 inter- 
ference. The amplitude for the process in eq. (2.6) with n — 1 is 

A = M» q2 _ ml \. mvTv (-gr + ||) WHW (Vv, - ^ 75 ) w(fe ) , 

(2.15) 
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where my and IV are the mass and the width of the vector boson respectively, and 
M M is the amplitude for the process 

a{P x ) + b(P 2 ) — V(q) + X(x) , (2.16) 

p being the Lorentz index associated with V; the polarization four- vector of V is not 
included in M M . From eq. (2.15) we get (neglecting lepton masses) 

12 _ + QnQu/ml) {-g pa + g p q a /m v ) 

\2 



x Tr[(Vyj + A 2 vl - 2V vl A va5 ) hYhl"] • (2-17) 



spins 



We now consider the narrow width approximation IV — > 0. We have 

1 7T 



5(q 2 -m 2 v ). (2.18) 



(g 2 - m 2 ,) 2 + (m y rV) 2 m v T v 
The 5 function, which puts the vector boson on shell, allows us to write 



+ ) = I>^r , (2.i9) 



where e\ are the polarization four- vectors of the vector boson. Using eq. (2.19), 
eq. (2.17) becomes 

5>! 2 = £ M\P\\'My 5 (q 2 - m 2 v ) , (2.20) 

spin ^ ^ AA' 

where we defined 

M A = M^, (2.21) 

which is the amplitude for the process of eq. (2.16) for a given vector boson polar- 
ization A. We also define 

Pxx , = F 2 Tr[(V 2 l + A 2 vl -2V vl A va5 )hfxhtx'} (2-22) 

which is, apart from the normalization, the decay density matrix 2 of the vector boson. 
This quantity can be explicitly computed; here, we only present it in the form of a 
diagonal matrix 

Paa' = {U P D U*) xy , (2.23) 

where 

p D = 2m 2 v F 2 diag (0, (V vl - A vl ) 2 , (V vl + A vl ) 2 ) . (2.24) 



2 The density matrix is usually denned as the transpose of that in eq. (2.22). See e.g. ref. [4]. 
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The cross section for the production of a lepton pair in the narrow width approxi- 
mation is therefore 

= h^k- E MA M, 6 ^ - "*) 

AA 

da 2 

x ^- d$i+i*(Pi + P 2 ; q, x) d$ 2 (g; fci, fc 2 ) . (2.25) 

Owing to the hermiticity properties of the density matrices, and to the explicit form 
of eq. (2.24), eq. (2.25) is a positive-definite quadratic form in the space of the spin 
indices of the vector boson. Thus 

ia '> £ i^r v t i*i s d * 2 • (2 - 26) 

where 

|M| 2 = ^ M A M A * = (MU) (mu)* = M"M " f-g^ + ^ j , (2.27) 

with g 2 = my. Eq. (2.26) cannot be used as an upper bound for the matrix element of 
the process (2.6), since the measures on the two sides are different. Using eq. (2.18), 
we can however easily reinstate the q 2 integration by inserting 



/ 



7i (q 2 - m v ) 2 + imyYy) 2 
on the r.h.s. of eq. (2.26). Furthermore, from eq. (2.24) we obtain: 

max (pf A ) = 2m 2 v F 2 (V vl + A vl ) 2 , (2.29) 

A 

which holds since V V iA V i > regardless of the identity of the lepton /. Therefore 

lyup, da fl 2m 2 v F 2 (V vl + A vl ) 2 \M\ 2 
2s ^ 1 d$ 2+v ~ (q 2 - m 2 v ) 2 + (m v T v ) 2 2s ' 1 ' ' 

which strictly speaking holds only when q 2 = my, since all results in this section are 
formally derived in the limit IV — > 0. More details on this, and the reason for keeping 
a formal dependence on q 2 in eq. (2.30), will be given in appendix A. Eq. (2.30) is 
the main result of this section. It states that, in the narrow width approximation, 
the lepton-pair cross section has an upper bound, which is a universal factor times 
the cross section for the production of the parent vector boson 

dav =-\M\ 2 . (2.31) 



2s 
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2.3 Top decay 

Here, we consider the decay of a top quark 

t{p) — > W + (q) + b(r) — > Z+(fci) +1/(^2) +6(0; (2-32) 

the treatment of the decay of an antitop is fully analogous. Other top quarks may 
be present in the final state, but their decays are of no interest for the moment, and 
will be ignored. The amplitude for the process in eq. (2.12) is 

A = u{r)V tb ^{l - 75) 2 **t?- r M 
2V2 p 2 -mf + im t Y t 

q 2 - m 2 w + tm w T w K V 2V2 K ' K ' 1 ' 

where M is the amplitude for the process 

a(Pi) + 6(P 2 ) — > *(p) + ^(x) , (2.34) 

except for a spinor Tt(p), which is not included. Therefore, M = Tu(K), with T a 
combination of 7 matrices, and the four-momentum of a fermion entering the hard 
scattering. By squaring eq. (2.33) we get 



\A 



4 IT/ 1 2 

2 <7 w rt&l 



64 (p 2 - m 2 ) 2 + (m t r t ) 2 (g 2 - m 2 w ) 2 + (m^IV) 2 
x w(r)7 M (l - 75) {j> + m t ) MM'- (/} + m t ) (1 + 75)7^) 

x u(fc 2 )7 M (l - TsMfciMfci) (1 + TsfrXfe) • (2.35) 

Following what was done in eq. (2.17), we now consider eq. (2.35) in the narrow 
width approximation Y t — > 0, i.e. we make the replacement 

t-o — 7 — wrz — ► S (p 2 - m 2 ) . (2.36) 

(p 2 - m 2 ) 2 + (m t T t ) 2 m t T t ^ tJ 1 ; 

Thanks to the on-shell condition introduced in this way, we can use the analogue of 
eq. (2.19) 

p' + m t = J2Mp)Mp), (2-37) 

A 

which in turn suggests introducing the quantity 

M x = u x (p)M M A * = M* 7 Vp), (2.38) 

which is the analogue of eq. (2.21), and is the amplitude for the process of eq. (2.34) 
for a given top polarization A. Upon summing over the spins of the final-state leptons 
and b quark, eq. (2.35) can be cast in the same form as eq. (2.20): 

E I A \ 2 = ^f t E MxPxyMy 5 {p 2 - m\) , (2.39) 

spin XX' 
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with 

4 I y 1 2 -j 

Pay = w 1 a tb y~ 2 TyTTl r — V2 ^K 1 ~~ 7s)#27^i7p] 

x %(p)(l + 7 5 )f?V«A(p). (2-40) 

This is the decay density matrix for the top quark, the analogue of eq. (2.22). We 
can now proceed exactly as was done in sect. 2.2, and therefore we must compute 
the decay density matrix, diagonalize it, and find the largest of the matrix elements 
so obtained. An explicit computation leads to 



n *<& \V, 



Using eq. (2.41) and reinstating the integral in dp 2 using the analogue of eq. (2.28), 
we finally arrive at 



< 



2s ^'"" d *™* ^2 _ m 2 w y + {m w T w )^({p^ - mlY + (m t r t ) 2 ) 2s 

(2.42) 

which strictly speaking holds only when p 2 = m 2 . Eq. (2.42) is the analogue of 
eq. (2.30), and expresses the upper bound on the matrix elements for the production 
of Ivb in terms of the matrix elements for the production of a top quark 

litrh^- (2 - 43) 

In contrast to eq. (2.30), the bound of eq. (2.42) is not a constant over the phase 
space of the particles emerging from top decay, because of its dependence on (r • A; 2 ) 
and (p-ki). This helps to increase the efficiency of event generation in the context of 
an unweighting procedure, but in order to avoid any biases the phase-space must be 
sampled in such a way as to reproduce exactly the q 2 -, (r-A^)-, and (p-fci)-dependences 
of the bound. An alternative approach is that of finding a constant larger than or 
equal to the bound, which can be done by finding the maximum of the combination 
of dot products 

D = (r-k 2 )(p-k 1 ) . (2.44) 

Using the top rest frame to perform the relevant computations, it is a matter of 
simple algebra to obtain 

D<D m M) = \i/ 1 ' 2m ir\' 2 (2.45) 
[q {mf - q )/A -± < q 2 < mf . 

Note that m^/16 > q 2 (m 2 — q 2 )/A in the whole q 2 range, and therefore one can always 
set -D max = m£/16; this is seen to lead to a very marginal degradation of unweighting 
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efficiency. We have therefore 

dowb < %4 |^&| 2 -Pmax(g 2 ) da t ^ ^ 



d$ 3+i* ^2 _ m 2 w y + ( mw r w ) 2 ) ((p 2 - m 2 ) 2 + (m^) 2 ) 
2.4 Multiple decays 

It is easy to generalize the formulae derived in the previous sections to the cases in 
which one is interested in the decay products of several vector bosons and/or top 
quarks. Consider for example the process of eq. (2.6). An equation identical to 
eq. (2.15) holds, with the formal replacements 

M M — ► M^ 1 -^ , (2.47) 

J- J. n 2 — m 2 . -I- im-,r Ft, ' 



q 2 - my + im v T v ^1 g 2 _ m 2^ + ^^p^ 
FyuikjY (V vl - A w7 s) u(fc 2 ) — II ^"(Afc-ifr" (V^,, - A Vihlb ) v{k 2i ) . (2.49) 



i=i 

The analogue of eq. (2.25) features the quantity 



£ . . . (MU! . . . U n ) PiAiA'j • • • Pn\ n \' n (MU, ...U n ) * , (2.50) 

\ / Ai...A„ V / A,...A„ 

AlAj A n A n 



which results from the simultaneous diagonalization of the spin density matrices of 
the n vector bosons. This allows one to use eq. (2.29), and to proceed as in the 
previous section. We therefore arrive at 

da hh...U n < (yr Zm 2 Vi F 2 z {V v ^ + A v ^) 2 \ da Vl ... Vn 
<«W ~ [} = \ (Qf - KJ 2 + (rnvT Vi ) 2 J d$ n+1 * ' l ' > 

where dav 1 ...v n is the cross section for the process of eq. (2.5), all the vector bosons 
being on-shell. 

Along the same lines, eq. (2.46) can immediately be generalized to the case of 
the decays of several top and antitop quarks: 

d(Tl 1 u 1 b 1 ...l n i/ n b n ^ 



4gi\V tb \ 2 D max (q 2 ) \ da h ... tri 



i ((tf - m 2 w ) 2 + {m w T w ) 2 ) {{p 2 - m 2 ) 2 + {m t T t ) 2 ) ) 



. (2.52) 



Obviously, eqs. (2.51) and (2.52) can be combined for the simultaneous presence of 
vector bosons and top quarks in the final state. 
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3. Angular correlations in MC@NLO 

As mentioned in the introduction, the straightforward way to predict correctly all 
features of lepton spectra is to include in the computation the leptonic matrix el- 
ements, for example as done in MC@NLO version 3.3 [5] for the cases of single-V^ 
or VH production, or in refs. [6, 7] for the case of top quark decay in parton-level 
pure NLO computations. We remind the reader that, in the context of MC@NLO, 
parton-level cross sections (which are obtained by suitably modifying those which 
enter pure-NLO computations) are first integrated over the phase space of the final- 
state particles. The information gathered in this integration step is then used in 
the event-generation step, whose aim is that of obtaining a set of kinematic config- 
urations (the hard events), which are subsequently showered by the parton shower 
Monte Carlo. We also point out that the same integration-and-generation structure 
is used by POWHEG (although the cross sections integrated in the two formalisms 
are not the same). The integration time increases rapidly with the number of final- 
state particles; there is a corresponding decrease in the efficiency of the generation of 
hard events. This is the reason why it is interesting to find alternative ways to pre- 
dict angular correlations in large-multiplicity processes. We stress that in principle, 
the implementation in MC@NLO (or POWHEG) of a process with correct angular 
correlations is identical to that of the same process without such correlations. The 
problem is a practical one, namely that production angular correlations require the 
knowledge of the lepton matrix elements, and the increased multiplicity with respect 
to the undecayed matrix elements entails loss of accuracy and generation efficiency. 
The strategy we propose in this paper starts with the following steps. 

1. Integrate the undecayed matrix elements. 

2. Generate hard events using the results of the previous step; thus, vectors bosons 
and/or top quarks will be present in the final state, but not their decay prod- 
ucts. 

3. For each hard event, generate (massless) lepton (and b quark, in the case of 
top decays) four-momenta, uniformly within the decay phase space(s) of the 
corresponding parent particle(s). 

4. Compute the lepton matrix element using the four-momenta obtained in step 
3, and the undecayed matrix element using the four-momenta obtained in step 
2. Generate a flat random number r. If the lepton matrix element, divided by 
its upper bound as given in eqs. (2.51) and (2.52), is smaller than r, throw the 
lepton four- momenta away, and return to step 3. 

5. Otherwise, replace the vector bosons and top quarks by the set of their decay 
products. The resulting kinematic configuration is the leptonic hard event that 
can be showered by the Monte Carlo. 
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Steps 3 to 5 constitute a standard hit-and-miss procedure, which guarantees that the 
lepton spectra reconstructed with the four-momenta of the leptonic hard event (and 
subsequent shower) will be identical to those computed by a direct integration of the 
leptonic matrix elements. 

It is clear that the integration step will be greatly simplified by this proce- 
dure: the number of phase-space variables relevant to the undecayed processes (2.5) 
and (2.10) is riu = 3{n+n x ) — 4, whereas riy = 3(2n+n x ) — 4 and n t = 3(3n+n x ) — 4 
for leptonic processes (2.6) and (2.12) respectively. On the other hand, one may 
doubt that the efficiency for producing leptonic hard events is larger than in the case 
of a straightforward integration of the leptonic matrix elements. In fact, the adap- 
tive integration performed in step 1 will only give information on the nu degrees of 
freedom of the undecayed processes. However, using the phase-space decompositions 
of eqs. (2.8) and (2.13), one associates the extra n v — riu = 3n and n t — riu = 6n 
degrees of freedom with the decay phase spaces. Since we are considering here only 
resonant diagrams, the leptonic matrix elements will be fairly smooth in these extra 
3n and 6n degrees of freedom, if the parametrizations of the decay phase spaces are 
properly chosen (the obvious choice of using the rest frame of the decaying particles 
is also an optimal choice from this point of view). Therefore, all of the complications 
due to the presence of several peaks in the matrix elements are dealt with in step 
2. The unweighting performed in step 4 does not require any sophisticated numer- 
ical approach (i.e., a preliminary adaptive integration is not necessary) in order to 
achieve a satisfactory efficiency. 

For the procedure as outlined above to work, it is crucial that the leptonic matrix 
elements can be bounded from above by the undecayed matrix elements. In the 
derivations of sect. 2 we have assumed that the density matrix is positive definite, 
which is the case, and that the matrix elements involved can be expressed as the 
modulus squared of an amplitude. This is certainly the case in the context of a 
tree-level computation, but it is not true for all the contributions to an NLO cross 
section. In particular, the interference between virtual and Born amplitudes is not 
positive-definite in general. The modified subtraction procedure [1] introduced in the 
MC@NLO formalism also implies the presence of a second quantity which is possibly 
not positive-definite, namely the difference between the real matrix elements and the 
MC subtraction terms. The presence of non-positive-definite contributions is what 
prevents one from including angular correlations exactly to NLO accuracy in the 
context of the decay chain approximation, as anticipated in sect. 1. 

Before proceeding, we remind the reader that there are two classes of MC@NLO 
hard events, defined according to their kinematics: § (H) events have the same num- 
ber of initial- and final-state particles as Born (real-emission) contributions. Thus, 
the number of final-state particles of H events is equal to that of § events, plus one. 
For example, in W + W~ production (eq. (2.5)) we have (n,n x ) = (2,0) for § events, 
and (n,n x ) = (2,1) for H events. In tt and single-top production (eq. (2.10)), we 
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have (n,nx) = (2,0) and (n,nx) = (1,1) for § events, and (n,nx) = (2,1) and 
(n, nx) = (1,2) for EI events respectively. POWHEG (and, for that matter, any 
NLO computation) also outputs § and H events. 

We now extend the procedure proposed in points 1 to 5 above to the case of 
NLO computations matched to parton shower simulations, as follows: 

• Steps 1 and 2 are unchanged. 

• For each § event, go through steps 3 to 5, using Born-level results to compute 
lepton matrix elements and their upper bounds. 

• For each EI event, compute a quantity Q(M) as explained below and generate a 
random number r'. If r'< (7(H), go through steps 3 to 5, using real-emission re- 
sults to compute lepton matrix elements and their upper bounds. If r'> Q(M), 
define an S-type event with the projection Pe^s(H), and proceed as explained 
for S events above. 

The definition of a map Vu->s is a necessary condition for the matching between 
NLO results and parton shower simulations: for more details see e.g. refs. [2, 8]. 
This implies that such a map need not be defined specifically for the purpose of 
including angular correlations into MC@NLO or POWHEG. The quantity Q is a 
largely arbitrary smooth and continuous function, that assumes values between 
and 1, and tends to (1) in the soft /collinear (hard-emission) regions. The role 
of Q is simply to avoid computing real-emission matrix elements in the phase-space 
regions where they diverge. In the context of MC@NLO, functions with the same 
behaviour as Q need be introduced in order to ensure local cancellation between real 
matrix elements and MC counterterms (see e.g. app. A. 5 of ref. [1] and app. B of 
ref. [2]), and one obvious choice is that of setting Q equal to one of these functions 
(or to a combination of them). 

It should be clear that the proposal made here accounts for angular correla- 
tions to LO accuracy close to the soft and collinear regions, since there (/(EI) ~ 0, 
and therefore H events are projected onto § events, for which we only consider the 
Born matrix elements in the hit-and-miss procedure 3 . On the other hand, in the 
hard emission region only real corrections contribute to the cross section, and thus 
angular correlations are included exactly to NLO accuracy 4 . Angular correlations re- 
sulting from an MC matched to an NLO computation and implementing the method 
proposed in this paper have therefore the same or a better accuracy than LO-based 
Monte Carlos. We also stress that angular correlations are actually fairly close to 
those computed exactly to NLO, for two reasons. First, NLO corrections to spin 
correlations are generally small. Second, although virtual corrections and subtracted 

3 We remind the reader that the full NLO undecayed matrix elements are used in steps 1 and 2. 
4 One should bear in mind that radiation from the decay products is not included here. 
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terms are not positive definite, their angular correlations arising from the contribu- 
tions (if any) that are proportional to the Born matrix elements can be included 
exactly in the computation following the method proposed here, since both sides of 
eqs. (2.51) and (2.52) then get multiplied by the same factor. 

4. Results 

The approach described in the previous section has been adopted to include pro- 
duction angular correlations in MC@NLO in the cases of W + W~ production (since 
version 3.1) and of ti and single-t production (since version 3.3). In this section we 
present sample results for ti and single-i production, at the LHC (pp collisions at 
y/S = 14 TeV) and at the Tevatron run II (pp collisions at v^S* = 1.96 TeV). All 
the predictions given in this section have been obtained by using the MRST2002 
default PDF set [9], and by setting m t = 175 GeV and T t = 1.7 GeV. In the case 
of single-t production, we also reconstruct the accompanying jets, by means of the 
fc T -clustering algorithm [10], with d cu t = 100 GeV 2 . We include in the clustering 
procedure all final-state stable hadrons and photons. For the sake of simplicity, we 
force 7r 's and all lowest-lying 6-flavoured states to be stable in HERWIG. The jets 
are ordered in transverse momentum. 




100 150 200 o 1 

p T (i + A0(i + r) 



Figure 1: Transverse momentum of the lepton pair (left pane), and difference in azimuthal 
angle between the leptons (right pane), in production at the LHC. HERWIG results have 
been rescaled (by 0.3 on the left, and by the factor K = o" NLO /<7 LO on the right). 

We begin by considering ti production. We have studied, at the Tevatron and 
at the LHC, single-inclusive p T and rapidity spectra of the t and i decay products, 
and the correlations in transverse momentum, A</>, and invariant mass of the bb, 
l + l~, bl~, bl + , bu, and bv pairs. We have found that angular correlations have an 
almost negligible impact. We present in fig. 1 the only two observables for which 
these correlations have a visible effect, albeit barely so for p T (l + l~). On the other 
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hand, angular correlations are an important ingredient for the correct prediction of 
A0(/ + Z~), as shown in the right pane of fig. 1. It is interesting that about 30% of the 
difference between the LO prediction without angular correlations (dashed histogram 
- HERWIG) and the NLO prediction with angular correlations (solid histogram, 
overlayed with open circles - MC@NLO) is due to beyond-LO corrections. 

It is possible to specifically design observables which would be trivial if angular 
correlations were neglected. Typically, such observables are angular variables con- 
structed with the decay products of the top quarks, and measured in the rest frames 
of the parent particle. We have considered the distributions in cos#i, cos6 , 2 , and 
cos0, as defined in ref. [11]; in particular, is the angle between the direction of 
flight of / + and the direction of flight of l~ . The directions of flight are defined in the 
t and i rest frames respectively (see ref. [11] for more details). Results for cos0 are 
presented in fig. 2 for the Tevatron (left pane) and the LHC (right pane). Just as 
for A0(/ + /~), beyond-LO contributions are not negligible, and they tend to deplete 
(at the Tevatron) or to enhance (at the LHC) the LO predictions for the cos asym- 
metry. This behaviour is also found in the pure-NLO, parton-level study of ref. [11]. 
We have verified that, by neglecting angular correlations, the cross section depends 
trivially on 9%, 82 and 4>. 
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Figure 2: Opening angle distributions, as defined in the text, for ti production at the 
Tevatron (left pane) and at the LHC (right pane). HERWIG results have been rescaled 
by the K factor. The corresponding curves obtained by neglecting angular correlations are 
flat, and are not shown in the figure. 



Finally, we examine distributions for single-top production at the Tevatron. Be- 
cause both production and decay occur through the left-handed charged current 
interaction, one expects stronger production angular correlations than in top quark 
pair production. Indeed, angular correlation effects are clearly visible in the single- 
inclusive spectra of the top decay products. As in the case of tt production, it is 
possible to study angular correlations more directly by choosing specific observables. 
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These observables always involve the definition of a spin basis that leads to nearly 
100% correlation between the direction of the charged lepton from top decay and 
another experimentally-definable, channel-dependent direction [12]. For both s- and 
t-channel processes the optimal spin quantization axis lies, in the top quark rest 
frame, along the down-type quark attached to the vertex connected via a VF-boson 
to the top quark producing vertex. At LO that corresponds for the s channel to the 
beam-direction, while for the t channel this is most often the direction of the light 
quark jet against which the top quark recoils. 
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Figure 3: Angular correlations in single-top processes at the Tevatron: cos 9 in single-t/t 
production (left pane) , and cos x m single-t production (right pane) . Histograms without 
symbols are the sums of s- and t-channel contributions. 



Accordingly, we present in the left pane of fig. 3 the distribution in the cosine 
of the angle 9, defined as the angle between the direction of flight of the lepton 
emerging from top decay, and the axis of the hardest jet which does not contain a 
stable 6-flavoured hadron; the angle is defined in the rest frame of the top quark. This 
distribution has been shown in ref. [13] at tree level, and in ref. [14] at NLO using 
MCFM [6]. We have applied similar cuts as those in ref. [13], namely we required 
the decay products of the top to have 

p T (b)> 20 GeV, \r](b)\<2, (4.1) 
p T (/)>10GeV, |r/(/)|<2.5, (4.2) 
p T {u) > 20 GeV. (4.3) 

We also require the hardest light jet to have transverse momentum larger than 20 
GeV, and \rj(j)\ < 2.5. In this way, we obtain A = —0.35, where 

„ a(-l< cos9 < -0.1) -o-(-0.1< cos0 < 0.8) , , 

A = — = = (4 4) 

o(-l < cos6 < -0.1) +o-(-0.1 < cosfl < 0.8) v ; 

As can be seen from fig. 3, this result is due to the contribution of the t-channel, 
the s-channel having a very small asymmetry. We remark that the asymmetry is 
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also compatible with zero if spin correlations are switched off. It is important to 
notice that our results follow the same pattern (and are actually close numerically) 
of those of refs. [13, 14]. Although we did not carry out a comprehensive study, this 
fact implies that not only is the cos 9 asymmetry fairly robust when including higher 
order corrections, but it is also stable when passing from a parton-level description 
such as that of refs. [13, 14] to a more realistic hadron-level description such as that 
of MC@NLO. 

We conclude by presenting in the right pane of fig. 3 the distribution in the cosine 
of the angle x, which is defined analogously to the angle 9, except for the fact that the 
reference direction is chosen to be that of the antiproton beam (at variance with the 
case of cos 9, we have limited ourselves here to considering t production, rather than 
t + t production). As expected [12], the dominant contribution to the asymmetry is 
due in this case to the s-channel. An extremely small non-zero asymmetry may also 
be visible in the case in which angular correlations are not included; we have verified 
that this is an artifact of the cuts adopted in the present analysis. 

5. Conclusions 

We have presented a method for the efficient inclusion of angular correlations due 
to production spin correlations in Monte Carlo event generators. The method has 
been demonstrated in detail for vector boson and top quark decays, but it is in fact 
quite general, relying only on the fact that the matrix elements do not contain sharp 
features that would lead to unacceptably low efficiency. The method is exact, and 
equivalent to what is currently implemented in LO-accurate event generators such 
as HERWIG. When the event generator is matched to NLO predictions, as is the 
case for MC@NLO and POWHEG, the resulting correlations are correct to LO in 
soft and collinear regions and to NLO elsewhere. The method has been implemented 
in MC@NLO for WW, ti and single-top hadroproduction and leptonic decay, and 
we have presented illustrative results for the latter two cases. These results show 
that significant correlations are present in suitably chosen observables. Version 3.3 
of MC@NLO implements off-shell effects only in the case of WW production. Future 
versions will include off-shell effects in top decay; also, vector bosons and top quarks 
decaying hadronically can be simulated using the formalism presented here, bearing 
in mind that NLO corrections to decays are neglected. 
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A. Upper bounds in vector boson production 

In this appendix, we present an alternative derivation of eq. (2.30). One introduces 
the quantity 

* ( v . %<t 



N u = M»- - g ; + ^L- , (A.l) 

with which eq. (2.17) becomes 

\A\ 2 = Fy Tr [(V^ + A vl - 2V vl A vllb ) MM*] ■ (A.2) 

spins 

Evaluating this in the rest-frame of the (virtual) vector boson, with the z-axis along 
the direction of the lepton 3-momentum, we find 

\ A \ 2 = 2( ? F v [( v vi + A vi) {^N 1 * + N 2 N 2 *) + AV vl A vl Im (N l N 2 *)] . (A.3) 

spins 

To establish an upper bound on this quantity, we note that 

2\\m{N 1 N 2 *) \ <2|A^| \N 2 \ < N X N U + N 2 N 2 * (A.4) 

and so 

\A\ 2 < 2q 2 F 2 (\V vl \ + \A vl \) 2 (A^ 1 A^ 1 * + N 2 N 2 *) 

spins 

= K^T^V^f* {M l M u + M 2 M 2 *) . (A.5) 
{q 2 -m 2 yf + {m v T v ) 2 K 1 

Note that M 1,2 in this expression are strictly off-mass-shell quantities: no on-shell 
approximations have been made at this stage. 

Now consider the production of a stable vector boson of mass my. Denoting the 
amplitude for this by A, we have (again in the vector boson rest frame) 

Y,\a\ 2 = m»m» (-g^ + w) 

spins V m V ' 1 2 =K 

= M '■.]!'■■ + M 2 M 2 * + M 3 M 3 * (A.6) 
where M M denotes the on-mass-shell value of M M . Therefore, as long as 

M 1 M 1 * + M-M-- < M X M U + M 2 M 2 * + M 3 M 3 * (A.7) 
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we have 



and hence 



E w . < ^n^ML E K (A . 8) 



(g 2 - m^) 2 + (m y r y ) 

spins spins 



da u ^ 2q 2 F 2 (\V vl \ + \A vl \) 2 day 



d$ 2 +i* (g 2 - m 2 /) 2 + (m v T v ) 2 cteV 

which is identical to eq. (2.30), given the fact that VyiAyi > 0, and that both 
equations are valid on-shell. 

Clearly, one may check whether the bounds given in eqs. (2.30) and (A. 9) are 
not violated in the case of off-shell vector bosons. This is indeed the case, provided 
that the off-shellness is not too large or too small (typically, this happens within 
±30IV of the pole mass). A good strategy is that of using eq. (2.30) for q 2 < rriy, 
and eq. (A. 9) for q 2 > my. However, one should bear in mind that in the case of 
off-shell particles the values of Bjorken x's, and hence of the PDFs, may change, thus 
potentially affecting the bound. 
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